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Rough dependence upon initial data exemplified by explicit 
solutions and the effect of viscosity 

Y. Charles Li 


Abstract. In this article, we present some explicit solutions showing rough 
dependence upon initial data. We also studies viscous effects. An extension of 
a theorem of Cauchy to the viscous case is also presented. 


Contents 


1. Introduction 

2. Explicit solutions with rough dependence on initial data 

3. Extension of a theorem of Cauchy to the viscous case 

4. Conclusion 
References 


1 

2 

8 

9 

9 


1. Introduction 

It is well-known that both Navier-Stokes and Euler equations are locally well- 
posed in the Sobolev space H s (s > | + 1) where d is the spatial dimension [4] |5j. 
In two dimensions (d = 2), the well-posedness is also global. The solutions are in 
the space C°([0, T), H 3 ), where T is either finite or infinite. The solution operator 
F l is a one-parameter family of maps on H 3 . For each fixed t, F* maps the initial 
condition u(0) to the solution’s value at time t, u(t). F t is continuous in zt(0). In 
[X] , a simple example was presented showing that the solution operator F 4 is not 
uniformly continuous in u(0) under the Euler dynamics. In |2j [3], the nowhere 
differentiability of F t under the Euler dynamics was proved. In [7], we derived the 
upper bound on the derivative of the solution operator under the Navier-Stokes 
dynamics: ||VE*|| < e cr 'TReVt+a 1 t w ]. lere a anc j ai are constants depending only 
on the norm of the base solution. In this article, we present explicit solutions 
to both 2D Euler and 2D Navier-Stokes equations. These solutions show that 
under the Euler dynamics, ||V.F t || = oo, and under the Navier-Stokes dynamics, 
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||V-F' t || —> oo as the Reynolds number approaches infinity. We name the above 
nature of dependence of the solution operator on initial data as rough dependence 
on initial data. 


2. Explicit solutions with rough dependence on initial data 


In this section, we present explicit solutions to the 2D Euler equations and 2D 
Navier-Stokes equations, which show rough dependence on initial data. Consider 
the following 2D Euler equations 

(2.1) dtu + u ■ Vu = — Vp, V ■ u = 0, 


and 2D Navier-Stokes equations 

(2.2) dtu + u • Vu = — Vp H-Au, V • u = 0, 

Re 

under periodic boundary condition with period domain [0,27r] x [0,27t], where 
u = ( 111 , 112 ) is the velocity, p is pressure, and the spatial coordinate is denoted by 
x = (xi,x 2 ). First we study the following simple explicit solution to the 2D Euler 
equations, 


(2.3) 


u 1 


E “ 3+7 sin[n(a ;2 - at)], u 2 = a, 

n =1 


where i < 7 < 1, and er is a real parameter. We view this solution as a solution 
in the space C°([0, 00 ), H 3 ), where H 3 is the Sobolev space. We select the Sobolev 
space H 3 because both the 2D Euler equations and the 2D Navier-Stokes equations 
are well-posed in H 3 . The same construction in this article can be easily extended 
to any Sobolev space H s (s £ 1). Denote by F t the solution operator of either the 
2D Euler equations or the 2D Navier-Stokes equations. For any fixed t , one can 
view F 1 as a map in H 3 , 

F* : u( 0 ) —> u(t). 

The initial condition of the solution (12.31) is 


Ul 


E 


1 

n 3+l 


sin(na;2), 


tt2 = cr . 


n= 1 

By varying a, we get a variation direction of the initial condition, 


du i( 0 ) = 0 , du 2 (0) = da, 

which leads to the directional derivative d^F* of F l : 

du\ — t r . ,, du 2 

= E ^ C0S ^ X2 - )]» iv = L 

n=1 


da 


The norm of d a F * is given by 


00 

WaP t \\ 2 H 3 = ^ + 27T 2 t 2 E(! + H 2 + ” 4 + n6 )^4+ 


4+27 


= 4t r 2 + 27T 2 t 2 


n= 1 


n =1 

1 1 


n 


4+27 n 2+27 


d 2^ 71 


2-27 


= 00 , (t > 0 ), 
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where the last series is divergent when i < 7 < 1. Thus the directional derivative 
d a F 4 does not exist. Therefore the derivative VF 4 does not exist in view of the 
fact that the norm of the derivative VF 4 is greater than or equal to the norm 
of the directional derivative d a F t . In fact, the solution operator F 4 is nowhere 
differentiable [2,. 

The corresponding solution of (12.31) to the 2D Navier-Stokes equations (12.21) is 

00 1 „ 2 

(2.4) ui = ^2 -g^e - ^ 4 sm[n(x 2 - crt)}, u 2 = a, 

n= 1 

which is also a solution in the space C°([0, 00 ), H 3 ). The directional derivative 
of F 4 is 


du\ 
~d.a 


E —t -n2-t r / ,\i du 2 

-^e cos [n(x 2 -at)\, — = 1 . 


The norm of d a F 4 is given by 


RF 4 ||^3 


= 4 - 7 T 2 + 27r 2 f 2 (1 + n 2 + n 4 + n 6 ) 

n =1 


IT 2 nr. 2 

-e -2- ^ 

n 6+2 7 c 


Let 

(2.5) 5 (£) = t 2 £e“^. 


The maximum of <?(£) is given by 


,/(£)= t 2 e-S« (l-|^) =Q, 


that is, 

( 2 . 6 ) 

where the maximal value of g is 


£ = 


i?e 

~2t’ 


t Re 1 
9 = —^~ e ■ 


Thus 


RF 4 || H 3 < (l + VR~eVi^j\ 


U(Uj|| ff 3. 


Notice that for the full derivative VF 4 , the upper bound is given as [7j, 


|VF 4 || ff 3^ ff3 < 


where er and a\ are two constants depending on the F 3 norm of the base solution. 
From (EU), let 



, the integer part of 


I Re 

~2t)' 


then 


1 Re 


I Re 


Re 


2 V 2 1 


< n < \ ——, when \ —— > 1 . 


2 1 


2 1 
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We have 


\\d*F*r H 3>4tH + 

27 T 2 t 2 (l + n 2 + n 4 + n 6 )-^e- 2 ^ 4 

> 47t 2 + 27r 2 i 2 n 2_27 e _2 ^* 

2-27 


> 47t + 27t i 


2.2 


1 Re 

2 V 2 ^ 


-1 


= ( 2^) 2 + 


{ 2^2 J 


I-7' 


> 


27r+ ^L 7 ^^ 


1-7' 


|| 5 cr -F' t ||ff 3 > v^ 2 tt + -^f 7 


y 2V2 y 

/ y/ty/Re\ 


l 2 v / 2 y 


II^<t^ ? IIh 3 —> 00, 


Thus 

As i?e —» 00, 
thus 

||VF*||jja -> 00. 

From another perspective, one can also obtain a lower bound for ||<9 (T F t || // 3. Fix 
a t > 0, for each n, the maximum (EH) specifies a Reynolds number Re^ (with 
t = n 2 ), 

Re {n) = 2 tn 2 , 

where 

g(n 2 ) — e _1 t 2 n 2 . 

Since 


RiHlffa >4 tt 2 


n 


2^ t (l+n +n +n 
>47r 2 + 27rW- 27 e- 2 ^ t , 


-- 2 / 


we have 


> 2 n\ 1 + — n 2-27 , when i?e = Re^ n \ 

V 2 e 

Therefore, as -Re* 71 ) —► 00 (n —>• 00), 

II^IIh- 00. 

More general solutions with the same property of rough dependence on initial 
data can be derived as follows. We expand the vorticity into a Fourier series 

w = ^ oj k e ik ' x 

ke z 2 /{o} 
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where 

u- k =TJk, k = (ki,k 2 ), x=(xi,x 2 ). 

When the velocity has zero spatial mean, the 2D Euler equations can be re-written 
as 

where 
(2.7) 

When p || q, 


Uk = ^2 q)uj p ojq 

k—p+q 


A(p,q) = -(\q\- 2 -\p\- 2 ) 


pi qi 
P 2 92 


A(p , q) = 0. 

This implies that for any k £ Z 2 /{0}, 


( 2 . 8 ) 


UJ 


= E 

n£Z/{ 0 } 


0 in(k-x ) 


is a steady solution to the 2D Euler equation (in vorticity form) | 6 j, where C n ’s are 
complex constants, and C_„ = C n . In terms of velocity, this solution has the form 


“i = E 

n£Z/{ 0 } 


in(k-x) 


n\k\ 2 


u 2 


_ _ ikiC n J, n (k-x) 


nGZ/{0} 


n|A :| 2 


Under a translation in velocity, this solution is transformed into the following time- 
dependent solution, 

‘2^ti „in\k 1 (xi — ait)+k2(x2—<T2t)] 


n\k\ 2 


Ul = <Ji + ^2 

n£Z/{ 0} 

ll 2 = <72 — V 

n\k\ 2 


(2.9) 

n£Z/{ 0} 

where G\ and a 2 are real constants. By choosing 


( 2 . 10 ) 




we get a solution which has the same property of rough dependence on initial data 
as the solution (12.31) . The initial condition of the solution is 


ui ( 0 ) = cti + ^2 


ih 


n£Z/{ 0} 


n s +7\k\ 


2 _ in(k\x\-\-k2X2) 

C i 


u 2 


( 0 ) = cr 2 - E 


ik\ 


r 3 + 7 |fc | 2 


p in(kixi+k 2 x 2 ) 


n€Z/{ 0 } 

By varying ay (or a 2 ), we get a variation direction 

du\(0) = dai, du 2 (0) = 0 , 











6 


Y. CHARLES LI 


and the corresponding directional derivative d Gl F t of the solution operator F* is 


du\ 
da± 

du2 

hr 

The norm of d ai F t is given by 


_ 1 I \ ' t ^1^2 m[fci(xi— <rit)+k 2 (x 2 — o- 2 t)l 

dax ^ n 2+7 \k\ 2 

ne Z/{0} 1 1 


9(71 ^ n 2 +^ \k\ 2 

nez/{o} 1 1 


i h 2 

L "'1 „in\ki(xi — p-it~)+fc 2 (x 2 —<72t)l 


fc? 


11 -F* 11 //3 = 4tt 2 + 47r 2 t 2 |^2 5Z 


l + n 2 |fc| 2 + n 4 |fc| 4 + n 6 |/c| 6 


n£Z/{0} 

1 

n 4 + 2 7 “ 


when f > 0 and k\ ^ 0. 

Under viscous effect, the corresponding solution of (12.91) to the 2D Navier-Stokes 
equations m is given by 


\ ' lk2C n n 2 \k 

u\ = a\ + y. —mire Re ”e 


2 °n -- " t in[k 1 (x 1 -<j 1 t)+k 2 (x 2 -<J2t)] 


n^L/ { 0 } 


n|fc | 2 


( 2 . 11 ) 


u 2 


= (J 2 _ V' ikiCn c _ " J,); 1 t c in[fc 1 (a 1 -er 1 t)+fc 2 (a: 2 -<T 2 *)l 

nlfcl 2 

ne Z/{0} 1 1 


We still choose C n as in (12.101) . The directional derivative d ai F t is given by 

Ul = 1 i \ ' t k\k2 n IfJ t i n lki(xi-crit) + k2(x2-cr2t)] 

da! ^ n 2+7 \k\ 2 

nez/{0} 1 1 


t fc? ™ 2 I^I 

-L 


dlt2 

9(Ji n 2+7 |/c| 2 

nez/{o} 1 1 

The norm of d (Jl F t is given by 


in[k 1 (x 1 -a 1 t)+k 2 (x2—T2t)] 


=47r 2 +47r 2 t 2 y^ ( 


1 + n 2 |&| 2 + n 4 |fc| 4 + n 6 |fc| 6 


n£Z/{ 0} 

2 -w 


n 


e 2 Re 


n 6+2 7 

By the result on the function g(£) (12.51) with £ = n 2 |fc| 2 , we get 


R^IIhs < 

From (12.61) . let 




W(Uj|| ff 3. 


n = 


1 Re 

\k\ V 2 f 


, the integer part of 


1 fRe\ 
Jk\ \ 2t)’ 
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We have 


\\da 1 F t \\jj 3 > 47T + 

(! + ^l^l 2 + n4 | fc | 4 + - 6 N 6 )^ 

> 47t 2 + 47r 2 t 2 fc 2 |fc| 4 n 2_27 e -2 «« * 




> 47 T 2 + 47 T 2 t 2 fc 2 |fc| 4 


1 1 Re 

2 ]k\\l ~2t 


2-27 


-1 


= (2tt) 2 + 


^fci|fc| 1+7 t 7 (VtVRe\ 

\/e 


1-7' 


1 

> - 

“ 2 


Thus 

As -Re —► 00, 
thus 


O 7 r 

27 T+^=fc 1 |fc| 1 + 7 t 7 

Ve 


V^TI 


\ 2 V 2 y 

/ y/t\/Re\ 


1 - 7 ' 


R.Fitfa >^277 + ^—^rki\k\ 1+1 f 1 

Ve 


-» 00, 


IVR‘11^3 —»• 00. 


/ \fi\fRe\ 

l^TT ) 


1-7 


Now we go back to the formula m, when \p\ = |g|, 

Mp, q) = 0 . 

This implies that 

w = ^ Cke lk ' x , (c is a constant) 

|fc|=c 

is a steady solution of the 2 D Euler equations (in vorticity form) [6], where Ck’s 
are complex constants, and C-k = Ck■ In terms of velocity, this solution has the 
form 

"i = E 

|fc|=c 1 1 


U2 


= -£ 


|fc|=c 


RcCfc ifc , 

l*f 


Under a translation in velocity, this solution is transformed into the following time- 
dependent solution, 

mi = en + V 

7/0 — /To — V i[fcl(zi-<Tit)+fc 2 (a; 2 -g 2 «)l 

2- 2 Ifcl 2 

|fe|=c 1 1 


( 2 . 12 ) 
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where a\ and <72 are real constants. A special case of (12.1211 was given in [lj. As 
shown in [l], by choosing 01 (and/or < 72 ) to be -k k\ (and/or k-fl) to be n, and C'k to 
be n 1+s (s £ M), we end up with a sequence of solutions. As n —> 00 , this sequence 
of solutions shows that the solution operator is not uniformly continuous in H s . 

3. Extension of a theorem of Cauchy to the viscous case 

Let £(x, t) be the fluid particle trajectory starting from x at the initial time 
t = 0 , 

(3.1) d t £( x , t ) = u(£(x, t),t), £{x, 0) = x. 

For any fixed t, one can think £(x, t) as a map of the fluid domain. By the Liouville 
theorem 

det (V x £(z. t)) = det (V x £(®, 0)) e tv f u = 1, 

since 

• u = 0, det {V x f;(x, 0)) = 1. 

That is, £ is a volume-preserving map in time, and any volume in the fluid domain 
is preserved during the fluid motion. The following is a theorem proved by Cauchy 

mm- 

Theorem 3.1. The vorticity along a fluid particle trajectory under the Euler 
dynamics, evolves according to 

(3.2) u}i(£(x,i),t) = &j(x,t)wj(x,0). 

The equation (13.21) has the same form with the variation equation 

S&(x,t) = £i,j(x, t)Sxj. 

It is not easy to extend the equation (13.21) to the viscous case, while an alternative 
form can be extended to the viscous case. One can introduce the tensor version of 
vorticity 

-li'j - C., ,'y Uj , . 

Then the theorem of Cauchy takes the following form [2] 

Theorem 3.2. The matrix vorticity along a fluid particle trajectory under the 
Euler dynamics, evolves according to 

(3.3) t)^kj — ^lij{x, 0). 

Equation m is a consequence of the equation 

At (£m,iHmfc(£k)£fc,l) = 0. 

The classical form of the Cauchy theorem (13.21) can be obtained from (13.31) when us¬ 
ing the u> variable. Equation (13.31) plays a key role in proving the non-differentiability 
of the solution operator for Euler equations [2|. Equation (13.31) can be extended to 
the viscous case. 

Theorem 3.3. The matrix vorticity along a fluid particle trajectory under the 
Navier-Stokes dynamics, evolves according to 

(£> ^')^k,j') = ~^^^m,i^^mk (£5 ty£k,j * 


( 3 . 4 ) 
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PROOF. Starting from the Navier-Stokes equations 

dt'U'i + Hj'U'iyj — P,i ”1” "75 ^^ii 
He 

we get 

(3.5) + Uk^Qkj + Uk^ij,k = . 

He 

From (ED, we have 

(3.6) dt£,i,j = Ui,k£k,j‘ 

The material derivative of Clij(£,t) is given by 

(3.7) d t to%j(€, t ) = d t Qij + u k ^ij,k- 

Now we can use m and ED to calculate the material derivative 

(£ra,i^rafc (£>*)£fc,j) = 6 ,1^771,l^mk£,k,j 

[Ot^mk + HZ^rafc,z]£/c,j "F £m,iQmk'U'k,l£,l,j 

— £,171,1^1,m^lk^kjj H - ^m,i [dt^mk + + £,m,i^mlUl t k£,k,j 

= T 4“ ^mL^l,k T ? 

by rearranging the dummy indices. Using (13.511 . we get (13.411 . □ 

4. Conclusion 

By constructing explicit solutions to 2D Euler and 2D Navier-Stokes equations, 
we showed the rough dependence of their solution operators upon initial data. The 
effect of viscosity is also studied. Finally we presented an extension of a theorem 
of Cauchy to the viscous case. 
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